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Abstract 

We study the effects of virtual leptoquarks on charged current and neutral 
current processes at the ep collider HERA. We present the areas of parameter 
space that can be excluded at HERA by searching for deviations from Standard 
Model expectations. The best results are obtained by examining the ratio of 
neutral current to charged current cross sections, R = ctnc/^CCi where, with 
200 pb _1 of integrated luminosity for unpolarized e~ and e + beams, HERA can 
search for leptoquarks with masses up to ~ 800 GeV, with leptoquark coupling 
strengths of order a em . 
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Many theories which go beyond the Standard Model (SM) are inspired by the 
symmetry between the quark and lepton generations and try to relate them at a more 
fundamental level. As a result, many of these models contain new particles, called 
leptoquarks, which naturally couple to a lepton-quark pair. Examples of such theories 
include[|l| models with quark-lepton substructure, the strong-coupling version of the 
SM, horizontal symmetries, grand unified theories based on the gauge groups SU(5), 
SO(10), as well as the Pati-Salam SU(4), and superstring inspired E 6 models. In all 
these theories, leptoquarks carry both baryon and lepton number, and are triplets 
under SU(3)c; their other quantum numbers (e.g., spin, weak isospin, and electric 
charge) vary between the different models. These particles need not be heavy; in fact, 
leptoquarks can have a mass < 100 GeV and still avoid[Q conflicts with rapid proton 
decay and dangerously large flavor changing neutral currents. This is particularly 
true in models where each generation of fermions has its own leptoquark(s) which 
couples only within that generation. Models where the leptoquarks can possess flavor 
non-diagonal couplings are much more constrained. 

By its nature, the high-energy ep collider HERA is especially well-suited to 
study leptoquarks. Direct production^ [| can occur via an s-channel resonance at 
enormous rates with distinctive peaks in the x-distribution. Indeed, HERA exper- 
imental searches are expected^ to reach a discovery limit (with a 5a signal and 
an integrated luminosity of 200 pb _1 ) of ~ 250 GeV when the leptoquark coupling 
strength is of order 0.01a em , and up to the kinematic limit if the coupling is equal 
to a em . If leptoquarks are too massive to be produced directly at HERA, perhaps 
they can be detected through their indirect effects via virtual exchange by searching 
for deviations from SM expectations for certain processes. This is similar in nature 
to the detection§] of the SM Z-boson at PEP/PETRA. Several authorsg, § have 
examined such effects on the neutral current asymmetries that can be formed with 
polarized electron beams, and have found that departures from the SM are small, 
even for leptoquarks of low mass (e.g., ~ 400 GeV) and large couplings. Here, we 
systematically investigate these possible indirect effects on charged current as well as 
neutral current processes. We find that the best results are obtained by examining 
the ratio R = unc/vcc, where discovery limits can reach leptoquark masses of order 
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800 GeV for electromagnetic coupling strengths (with 200 pb -1 of integrated lumi- 
nosity per e + , e~ beam). Taking the ratio of neutral to charged current cross sections 
is also advantageous because several systematic uncertainties, as well as those from a 
lack of detailed knowledge of the parton distributions, will cancel. Such a ratio has 
historically played a great role|| in our knowledge of SM interactions from traditional 
low-energy neutrino scattering. As such, it is only natural that such a ratio be used 
to probe new interactions in the high-energy regime at HERA. 

Present bounds from the LEP experiments |9|] on a leptoquark mass are m > 
Mz/2] single leptoquark production from Z decay could extend|l0| the LEP search 
reach to masses of ~ 70 GeV. It is also possible, using the high statistics available 
at LEP, to study the virtual effects of leptoquarks via flavor changing Z decays JIT]. 
Any such limits will, however, be more model dependent that those quoted above. 
At hadron colliders, leptoquarks may be produced [0 either singly through their 
unknown q£(LQ) couplings via qg fusion, or in pairs by gg fusion and qq annihilation. 
UA2 has placed [jn| limits on scalar leptoquark pair production by searching for the 
possible final states e + e~ + 2 jets and eu + 2 jets with the result that m > 67 GeV 
at 95% C.L. assuming the branching fraction B(LQ — > eq) = 50%. As one varies this 
branching fraction from 10% to 100%, the UA2 bounds range from the LEP result of 
m > Mz/2 up to m > 74 GeV. Constraints in the leptoquark mass-coupling plane 
may also be obtained by examining t-channel leptoquark exchange in the process 
e + e _ — > qq. Such effects are hidden at LEP I due to the Z-boson resonance, but 



data from PEP/PETRA allow [14] for leptoquark coupling strengths to be of order 
electromagnetic strength, or larger, for the mass range of interest here. 

For definiteness, we will concentrate on the leptoquark present in superstring- 
inspired Eq models^]. These leptoquarks, which we will denote by S, are scalar, 
charge —1/3, baryon number = +1/3, lepton number = +1, weak iso-singlets, and 
are the supersymmetric partner of the exotic color triplet fermion which appears in 
the 27 representation of Eq. Their interactions are governed by the Eq superpotential 
terms 

X L LS C Q + X R Su c e c + \'u c Sd c , (1) 
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where L and Q represent the left-handed lepton and quark doublets, respectively, and 
the superscript c denotes the charge conjugate states. The Yukawa couplings, i.e., 
the A's, are a priori unknown and for simplicity we set A^ = \ R in our numerical 
calculations. We make the assumption that the right-handed neutrino is too heavy 
to be produced at HERA and hence ignore any contributions from the A' term. If 
v c is light enough to be produced, the signature from its subsequent decay would be 
distinguishable from that of the SM charged current events considered here. The total 
leptoquark decay width (assuming the right-handed neutrino does not contribute) is 

r s = ^(2Ai + A^). (2) 

For calculational purposes, we parameterize the A's by 

Fl,r a, (3) 
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A L,R 
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and take F L — F R < 1. 

First we consider the contributions of leptoquark exchange on neutral current 
(NC) events. These e^q — > e^q and e^q — > e ± q events are generated by t-channel 
7 and Z exchange in the SM, and by s- and -u-channel leptoquark exchange. The 
Feynman diagrams responsible for these mechanisms are presented in Fig. la-b. The 
differential cross section for a left-handed polarized electron and an unpolarized pro- 
ton is 



da(e L p) 2iia 
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Eibl L (q) + b 2 LR (q)(l-y) 2 }x q (x) (4) 



dxdy sx 2 y 2 y q ^ d 

+ T,[bl R (q)+b 2 LL (q)(l-y) 2 ]xq(x) 



, s i {s — nil) xu(x) (l — y) 2 xu(x) 

+ F L b LL {u) sxyl 1 s > y > + 1 , U) 2 y 
{ l( s ~ m s) + (ms^s) 2 } [u - m 2 s \ 

1 , J t 2 xu(x) u 2 y 2 xu{x) Y\ 

+ 4 L{ L+ R > { [(s - m 2 ) 2 + (m s r s ) 2 ] + [u - m 2 } 2 f 
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and for a left-handed polarized positron and an unpolarized proton, 



da(e~Lp) 2na 
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dxdy sx 2 y 2 I ±J d 



Y.K L (<l) + b 2 RR (q)(l-y) 2 }x q (x) (5) 



q=u,d 

\ (1 — y) 2 xu(x) (s — mi) xu(x) 

+ F R b RR {u) sxy \ — f T — + — 2 —— 

{ [u - m z s \ [(a - m z s y + {m s Ts) 2 \ 

1 J? It? i E" 1 A f M V t 2 XM(x) \ N 

+ -.rR\r L + £ R ) < jt— + 



4 v I [« - ™IP [(a - ™l) 2 + (^r 5 ) 2 ] j / ' 

The corresponding cross sections for right-handed polarized e's are obtained by 

da{e%p) da(efp) 



dxdy dxdy 



with R. (6) 



Here q(x) and q(x) represent the quark and anti-quark distribution functions, s = xs, 
t = —Q 2 , u = —s + Q 2 , Q 2 = sxy, and x and y are the usual scaling variables. The 
functions bij are given by 

V vra (i - Mf) 2 + (M z r z )2 ' ^ 

with Q* being the electric charge of the incoming parton, C[ = sin 2 ^, ^ = 

— sin 2 and T[ L is the value of the fermion's third component of weak isospin. 
These expressions agree with those of Ref. H and also with the SM contributions of 



Ref. 16, 17 



The Feynman diagrams representing the SM and s-, u-channel leptoquark 
contributions to charged current (CC) events, e ± q — > uq and e ± q — > uq, are depicted 
in Fig. 2a-b. The differential cross section for a left-handed polarized electron is 
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+ F b sx f (f - m s) xu ( x ) ! (! ~ y? x d(x) 
L LL \ [(s - m 2 s ) 2 + (m s r s ) 2 ] [« - m|] 

Fft 2 x-u(x) ^ Flu 2 y 2 xd(x)\ 



4[(s - m| ) 2 + (m 5 r 5 ) 2 ] 4[u - m|] 2 / ' 
and for a right-handed polarized positron, 

,_, , I (1 — w) 2 xd(x) (s — m?) im(o;) | 

+ F L b LL sxy { + 1 



it — m 



SJ 



p-m|) 2 + (m 5 r s ) 2 ] 



Flu 2 y 2 xd(x) Fit 2 xu(x) 

77 o n n "I - 



4[w-m|] 2 4[(s-m|) 2 + (m 5 r 5 ) 2 ]; ' 

Here, the SM contributions agree with that in Ref. [|17]]. The leptoquark exchange 
diagrams also contribute to e R and e\ scattering, (note that there are no SM contri- 
butions to these processes), with 

da{e R j>) _ 2na 2 F L F R ( t 2 xu(x) | u 2 y 2 xd(x) \ 

dxdy sx 2 y 2 4 y [(i — m|) 2 + (msTs) 2 ] [u — m|] 2 J 

and 

da(e^p) 2-kol 2 F L F R ( u 2 y 2 xd(x) t 2 xu{x) \ 



dxdy sx 2 y 2 4 y [u — m|] 2 [(s — m|) 2 + (msTs) 2 ] 
For charged current processes, bn is defined as 



GpM 2 , t(t - Mjy) 



v^aTr (t - M 2 ,) 2 + (M^r 



\2 



(12) 



In order to ensure that the NC and CC events are cleanly separated and 
identified we impose cuts on the scaling variables x and y, as well as on the transverse 



6 



momentum of the out-going lepton. We have found that restricting x to lie in the 
range 0.1 < x < 1.0 satisfies [H]] the above condition, while yielding the greatest 
sensitivity to the indirect leptoquark effects. For a given value of x, we then restrict 
y to lie in the range 

max(0.1, ymin) <y< min(1.0, y max ) , (13) 

where, 

/ A l^cut\2 

(14) 



ymax,min 2 



xs 



In order to fully distinguish the two types of events, we impose the realistic cuts[p"9H 

p T (e)>5GeV (15) 

for the out-going e ± in the NC events, and 

^ T (i/) > 20 GeV (16) 

on the missing transverse momentum in CC events. In the circumstance that a more 
sophisticated detector triggering system is developed, we also ran a test case with a 
cut on the CC missing transverse momentum of foriy) > 10 GeV, and found that our 
results did not significantly change. We then integrate the above differential cross 
sections over these x and y ranges to obtain our results. 

We perform a x 2 analysis in searching for virtual leptoquark effects and com- 
paring them to SM expectations for various asymmetries, total cross sections, and the 
ratio R = <JncI°cc- We assume that most of the systematic uncertainties cancel for 
the asymmetries and in R, and only include the statistical errors for these quantities, 
which are given by 

6A=^£=-;i=, (17) 

for the asymmetries defined as A = (a a — cr / g)/(cr Q , + ap) with N a = Ca a being the 
number of events for a a for a value of the integrated luminosity C, and 

6R 11 



R VN^c VNcc 



where the '©' signifies that the errors are to be added in quadrature, and Nnc,cc is 
the number of neutral current and charged current events, respectively. In the case 
where we examine possible deviations in the total cross sections, we also include the 
systematic errors associated with the luminosity determination, giving 

5a 1 8C ,„„ N 

V = 7n 9 T' (19 > 

where we take|19[] the value 8C/C = 5%. Finally, we also use the sum of neutral 
current and charged current asymmetries in order to increase the statistics in searching 
for deviations from the SM. In this case, we calculate x 2 f° r the neutral current {Xnc) 
and for the charged current (xcc) processes separately, and then add 

Xtot = X 2 nc + X 2 cc ■ (2°) 
We then determine 90% and 95% C.L. deviations based on the value of Xtot- 



We present our results for the HMRS-B parton distributions of Ref. ||21|| . We 
have checked that the limits obtained from the asymmetries and the ratio of cross 
sections are insensitive to the choice of distribution functions as expected, but the 
bounds from the total cross sections can change by as much as 10% as the distribution 
function parameterizations are varied. In our numerical calculations, we use the 
valuesf2(J M z = 91.175 GeV, T z = 2.487 GeV, M w = 80.14 GeV, T w = 2.15 GeV, 
a{M z ) = 1/128, and x w = 0.233. 

We examine the effects using both unpolarized and polarized e beams. In 
the latter case, we define the polarized cross sections as 

a ± (+P) = i(l + P)a(etp) + \{1 - P)a(e±p) , 

&±{-P) =i(l-P)a(^) + l(l + PM^), (21) 

which takes into account that there will be only a finite degree of polarization. We 
take P = 80% and assume 125 pb _1 of integrated luminosity for each polarized 
beam. For the unpolarized case, we set P = in the above equations. Generally, we 
find that the use of polarized beams does not increase the leptoquark search reach 
over the unpolarized case (for a similar value of C). 
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We begin the discussion of our results by examining the six different neutral 
current asymmetries that can be formed from various combinations of the polarized 
e ± beams. Out of these six asymmetries, A£j£, which is defined as 

yields the best leptoquark search limits. These limits are shown in the leptoquark 
coupling-mass parameter plane in Fig. 3, represented by the dotted (dashed-dotted) 
curve in the upper left-hand corner which corresponds to the 90% (95%) C.L. Here, 
and in all of our results below, the discovery region lies to the left of the curves. 
It is clear from the figure that these types of measurements do not yield useful 
limits, reaching leptoquark masses of only mj~ 325 GeV, even with electromagnetic 
coupling strengths. This confirms the results of previous authors 0. [/]. 

We next discuss the limits that are obtainable from measurements of the total 
cross sections. Summing over both the neutral and charged current unpolarized cross 

sections, we find that the 90% (95%) C.L. discovery reach lies to the left of the 
solid (dashed) curves in Fig. 3, for 200 pb _1 of integrated luminosity per e~,e + 
beam. We find that the 5% uncertainty in the luminosity determination completely 
dominates the statistical errors, such that improved searches from this technique are 
not possible, even with substantial increases in the total integrated luminosity. Of 
course, measurements involving only cross sections are also severely hampered by 
the rest of the systematic errors as well as the uncertainties arising from the lack of 
knowledge of the parton distributions. 

A slightly more useful quantity in searching for indirect leptoquark effects is 
the (unpolarized) charge asymmetry A~ + , given by 

aye p) + aye + p) 

In Fig. 4 we show the 90% (solid curves) and 95% (dashed curves) C.L. discovery 
regions obtainable from the charge asymmetry for (a) 200 pb~* and (b) 500 pb" 1 of 
integrated luminosity per beam. Where possible, we present the results from neutral 
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current and charged current events separately, as well as those obtained from the 
sum of both types of events. As can be seen by comparing Figs. 3 and 4, the charge 
asymmetry can exclude a slightly larger region of parameter space than cross section 
measurements alone. 

It is, of course, possible to measure other quantities that are essentially inde- 
pendent of systematic uncertainties. We find two such quantities that provide much 
stronger search limits than any of the observables discussed so far. The first of these 
is the neutral current-charged current asymmetry, defined as 

A ZW = aNC - aCC , (24) 

&NC + °~CC 

and the second is the ratio R = unc/^cc- Here, a NC and a C c represent the total 
neutral current and charged current cross sections, summed over both electron and 
positron beams. The discovery regions in the leptoquark coupling-mass plane, based 
on measurements of these two observables, are displayed in Figs. 5 and 6 (note the 
scale change on the mj axis compared to Figs. 3 and 4). Figures 5 and 6a show 
the results from unpolarized e ± beams for integrated luminosities of 20, 200, and 
500 pb _1 per beam, while Fig. 6b presents the bounds obtainable from polarized 
beams, with 125 pb _1 per ef R beam. Comparing Figs. 6a and 6b, we see that 

although the total integrated luminosity (summed over all beams) is larger in the 
polarized case, it is the unpolarized measurements which yield the largest discovery 
region! Note also that the search regions resulting from measurements of R and A Z w 
are essentially equivalent. However, due to slightly different systematic uncertainties 
associated with identifying NC versus CC events, the asymmetry A zw may not be as 
clean a measurement as the ratio R. The region that can be explored with 200 pb _1 
per unpolarized beam is mj < 800 GeV for large leptoquark-electron-quark couplings 
(F ~ 1), and F > 0.13 for m s ~ 314 GeV, based on measurements of these two 
quantities. 

In summary, we have performed a careful analysis of the effects of virtual lep- 
toquarks at HERA. We found that measurements of the standard neutral current 
asymmetries using polarized beams, of the charge asymmetries, or of the neutral 
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current and charge current total cross sections are either plagued by systematic un- 
certainties or are simply not very sensitive to the parameters of the theory. However, 
it is possible, by examining ratios and asymmetries between NC and CC cross sec- 
tions, to eliminate the large systematic errors and to probe a sizable region of the 
leptoquark parameter space. In principle, HERA can search for leptoquarks up to a 
mass of 800 GeV for large leptoquark-electron-quark couplings, and for leptoquark- 
electron-quark couplings > 0.13a em for leptoquarks with mass mj > y/s. 
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FIGURE CAPTIONS 



Figure 1. Feynman diagrams contributing to the neutral current events (a) e ± q — > e ± q 
and (b) e ± q — > e ± q, in the Standard Model and with leptoquark exchange. 

Figure 2. Feynman diagrams contributing to the charged current events (a) e~q — > vq' 
and (b) e + q — > vq' in the Standard Model and with leptoquark exchange. 

Figure 3. Discovery region in the mg — F plane at HERA by neutral current plus charged 
current total cross section and by the asymmetry A^r, for the integrated lu- 
minosities as shown. The 90% (95%) C.L. discovery region corresponds to the 
area to the left of the solid (dashed) curves for the cross section measurements 
and to the dotted (dashed-dotted) for the asymmetry. 

Figure 4. Discovery region in the ms — F plane at HERA by the charge asymmetry, 
A~ + , for the charged current (Aqq), the neutral current (A^q) and the sum 
(Acc + A N + C ) for (a) C = 200 pb^ 1 and (b) 500 pb _1 . The area to the upper 
left of the solid (dashed) curves can be excluded at the 90% (95%) C.L. 

Figure 5. Search region in the m s — F plane at HERA by the neutral-charged current 
asymmetry, Azw for various values of C as shown. The area to the upper left 
of the solid (dashed) curves can be excluded at the 90% (95%) C.L. 

Figure 6. Discovery region in the mj — F plane at HERA from the ratio, R, with (a) 
unpolarized e ± and (b) polarized e ± beams, for different integrated luminosities 
as shown. The area to the upper left of the solid (dashed) curves can be excluded 
at the 90% (95%) C.L. 
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